Strangeness production from SPS to LHC 
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Abstract. 

Global strangeness production in relativistic heavy ion collisions at SPS and RHIC 
is reviewed. Special emphasis is put on the comparison with the statistical model and 
the canonical suppression mechanism. It is shown that recent RHIC data on strange 
particle production as a function of centrality can be explained by a superposition of 
a fully equilibrated hadron gas and particle emission from single independent nucleon- 
nucleoli collisions in the outer corona. 
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1. Introduction 



The enhancement of relative (to u,d quarks) strange quark production in high energy 



heavy ion collisions with respect to elementary collisions has been predicted long time 
ago to be a signature of the Quark Gluon Plasma (QGP) formation pQ. The idea 
was that chiral symmetry restoration favours strange quark production because of the 
reduced mass compared to its zero temperature constituent value. This abundant 
strangeness production could be observed provided that it survives hadronization, i.e. 
if the early produced strange quarks coalesce into hadrons without reannihilating. A 
specific prediction of such a mechanism is the enhancement of multiply strange particles, 
especially hyperons. 

These phenomena have been indeed observed: the ratio of newly produced strange 
to u,d quarks (the so-called Wroblewski ratio = (ss)/2((uu) + (dd)) shows about 
a factor 2 increase going from elementary to heavy ion collisions, (see fig. [1]) as first 
observed in ref. [2], and the hyperons shows a clear hyerarchical enhancement in central 
Pb-Pb collisions with respect to peripheral Pb-Pb and p-Pb collisions at top SPS energy 
(V$nn = 17.2 GeV), as observed by the WA97-NA57 collaboration [3J. Also, it seems 
that this ratio increases quickly in heavy ion collisions as a function of centre-of-mass 
energy going from 1 to few GeVs and stay constant thereafter. 

The big question is what is the origin of this observed strangeness enhancement. 
Is the original prediction of generation in the plasma and subsequent coalescence still 
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Figure 1. Wroblewski ratio in elementary and heavy ion collisions as estimated from 
statistical model fits. The superimposed dashed line is the predicted value from a fully 
equilibrated hadron gas, the solid line the interpolation from hadron gas with extra 
strangeness suppression 75. 



viable? Or, rather, are the excess strange quarks produced essentially at hadronization? 
Or, finally, is strangeness produced during an intense hadronic re-scattering stage, 
according to transport models ansatzl Before trying to answer these questions, it is 
necessary to address a preliminary very important issue, i.e. whether we have produced 
a completely equilibrated hadron gas or not. If we have a completely equilibrated hadron 
gas, strangeness content is completely determined and gives information on freeze-out 
state, but it is not a probe of earlier stage of the process. Solving this problem may 
have a considerable impact on our understanding of strangeness production in relativistic 
heavy ion collisions. 

2. Statistical model and strangeness undersaturation 

The main tool to probe the formation of an equilibrated hadron gas are the fits of 
the measured particle multiplicities or ratios to the statistical model, that is the ideal 
hadron-resonance gas. Many authors have performed such analyses trying to pinpoint 
the thermodynamical parameters of the hadron emitting source at the chemical freeze- 
out and their conclusions are vastly different in this respect. Some [I] conclude that a 
completely equilibrated hadron gas has been produced throughout the examined centre- 
of-mass energy range (from low AGS to RHIC), others [5] that this never occurs. The 
reason of such a dramatic difference in physical conclusions resides on one hand on data 
selection and, on the other hand, on parameter choice in fitting procedure. 
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Figure 2. Left panel: js a s a function of centre-of-mass energy in central heavy ion 
collisions (from ref. 0). Right panel: 75 as a function of centrality in Au-Au collisions 
at \/s NN — 200 GeV in central heavy ion collisions [7] . 



The conclusion that a completely equilibrated hadron gas is found relies on the 
use of midrapidity densities as experimental input compared with integrated yields of 
a single fireball at full hadrochemical equilibrium as theoretical model. The underlying 
idea is that, being the QGP fireball expected at midrapidity, its properties can be probed 
by using midrapidity densities. However, such an approach requires the existence of a 
region around midrapidity (a plateau) where the thermodynamical parameters do not 
vary much. Conversely, if the observed particle rapidity distributions are not sufficiently 
wider than that of a single fireball at the thermal freeze-out temperature, the use of 
midrapidity densities instead of integrated multiplicities artificially enhances heavier 
particles which have, in general, narrower rapidity width than lighter particles. This 
has two biasing effects: increasing the estimated temperature and enhancing the yield 
of particles carrying strange quarks which are generally heavier than non-strange ones 
[H] so that the strangeness undersaturation parameter 75 turns out to be approximately 
1 in these fits and essentially unnecessary. 

In fact, up to SPS energies, the rapidity distributions are not wide enough to 
allow the use of midrapidity densities. For instance, the pion rapidity distribution at 
\/s NN = 17.2 GeV has a dispersion width of 1.3, while that of a single fireball at 
rest with T = 125 MeV (the thermal freeze-out temperature) is about 0.8, hence not 
much smaller. Conversely, the measured width at RHIC at ^/s NN = 200 GeV is 2.1, 
which is reasonably larger than 0.8. Thus, the use of midrapidity densities allows to 
determine the thermodynamical parameters of the average fireball at midrapidity at 
RHIC energies onwards (roughly from 100 GeV), but not at SPS and lower energies. 
Therein, fits to full phase space yields provide a more appropriate, though amendable, 
estimate of the chemical freeze-out parameters. For recent studies including rapidity- 
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dependent chemical potentials see refs. PI [TUj. 

As a consequence, a strangeness undersaturation parameter 75 < 1 is needed to 
describe particle multiplicities in central heavy ion collisions. This parameter shows an 
increasing trend from AGS to RHIC, where it attains its maximal value 1 (see fig. [2]). 
Moreover, a 75 < 1 is also needed in peripheral collisions at RHIC for midrapidity 
densities, as shown in fig. EJ 

3. Canonical suppression 

It has been argued [UJ that the observed strangeness undersaturation (i.e. 75 < 1) at 
energies lower than SPS is owing to the so-called canonical suppression effect. Namely, 
strange particles are further suppressed in pp collisions and peripheral relativistic heavy 
ion collisions with respect to their expected yield in a grand-canonical ensemble (or 
thermodynamic limit) because strangeness is exactly vanishing within a small volume, 
called strangeness correlation volume (SCV), not necessarily coinciding with the global 
volume. Therefore, going from pp collisions to central heavy ion collisions through 
peripheral ones, one expects to observe a relative enhancement of strange particles 
due to approaching the thermodynamic limit, which is hyerarchical: Q yield increases 
faster than H which increases faster than A's or kaons. Yet, although this hyerarchy 
of enhancements is observed (see fig. [5]), neither SPS nor RHIC have observed the 
saturation which should be there if the SCV attains a sufficiently large value. In fact, 
this means that the SCV only reaches its saturation value (the one sufficient for the 
system to be essentially grand-canonical) at RHIC precisely in central collisions, where 
7s ~ 1. This would be quite a striking coincidence. Therefore, we think that canonical 
suppression is quite an unnatural explanation of the data, as already pointed out in 
ref. P2]. 

The best probe to investigate the phenomenon of strangeness undersaturation is 
indeed the meson. This is not an open strange particle, thus it is not canonically 
suppressed, yet, being a ss state, it must be 7I suppressed. Furthermore, meson 
has almost no feeding from heavier light-flavoured species and its production is entirely 
direct. 

It was pointed out quite early [13] that a statistical model with canonical 
suppression mechanism, i.e. with SCV as additional parameter, would have not been 
able to explain the deviation of the <fi meson yield from its grand-canonical value and 
this has been demonstrated in fits to NA49 multiplicities [H]. Recently, the STAR 
collaboration has measured the midrapidity densities of <fr meson very accurately and 
the observed pattern as a function of centrality clearly shows (see fig. [5]) that these 
do not scale linearly with the number of participants, rather the ratio to pp value 
increases rapidly at very peripheral collisions slowly saturating thereafter. This non- 
linear increase cannot be attributed to a variation of the chemical freeze-out temperature 
because this is astonishingly constant as a function of N p as shown in fig. [3] and proves 
that a genuine extra suppression related to the strange quark is needed, as also reflected 
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Figure 3. Chemical freeze-out as a function of centrality in Au-Au collisions at 
yfi NN = 200 GeV 0. 



in the 7s fitted value (see fig. [2]). 



4. Core-corona superposition 

Can we explain 75 in relativistic heavy ion collisions in more fundamantal terms? Some 
years ago [5] R. Stock proposed that 75 < 1 in global fits could be the effect of 
superposing a completely equilibrated hadron gas (75 = 1) originated from the core 
of the nuclear collision (i.e. the hadronization of the plasma) to a corona of single 
NN collisions where particle readily escape the interaction region. Since strangeness is 
largely suppressed in NN collisions with respect to the grand-canonical value while the 
temperature is almost the same as we know from pp statistical model analysis [HI [5], 
if the number of such single NN collisions accounts for a significant fraction of total 
particle production, a global fit to one hadron-resonance gas would actually find 75 
significantly less than 1. Indeed, this idea proved to be able to satisfactorily reproduce 
particle multiplicities in central C-C, Si-Si and Pb-Pb collisions at top SPS energy. 

This core-corona superposition mechanism has been invoked by Bozek few years 
ago [15J to reproduce the K/tt ratio as a function of centrality in Au-Au collisions and 
has been recently advocated in a paper by Werner [16] to be capable of explaining 
some more otherwise "mysterious" effects. A sharp superposition of a completely 
equilibrated hadron gas with NN collisions is indeed a zero-order approximation as the 
actual process is certainly more complex with those two extremes continuously linked 
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Figure 4. Number of single NN collisions divided by the number of wounded nucleons 
as a function of number of wounded nucleons according to Glauber Monte-Carlo. The 
solid line is the interpolation (TJ|. 



through intermediate steps and indeed in ref. [16j a more general concept of corona has 
been used, defined as a "dilute" peripheral region distinguished from the "dense" region 
in the core. Yet, this simple superposition scheme can be a very useful one to understand 
the physics of particle production. Accordingly, the rapidity density at midrapidity of 
any particle species is given by: 

= ^+w-««<|u. (i) 

where N s is the mean number of single NN collisions, Vq is the initial volume of the initial 
nuclear overlapping region, 5Vq is its thin outer shell where these single NN collisions 
occur, / is the growth factor (i.e. how much this volume expands up to chemical 
freeze-out) and dp/dy is the particle density per unit rapidity in the core relevant to 
a completely equilibrated hadron gas, i.e. with 75 = 1. Dividing by the number of 
wounded nucleons N w and the rapidity density in pp, we obtain a simple expression 
from ([I]): 



A+-?-(l-2A) (2) 



w 



where A is an unknown constant. This expression fulfills the constraint that both the 
left and right hand side ought to be 1/2 when N\y = 2 and N s = 1. Remarkably, for 

I In this work wounded and participant nucleons are synonymous 
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the (j) meson the constant A is independent of N\y because T is in fact independent of 
centrality (see fig. [3]) and does not suffer possible canonical suppression. Indeed, A is 
the asymptotic value of the normalized yield, when the number of participants becomes 
very large; hence, it is normally larger than 1/2 and the second term of the right hand 
side of eq. (J2J) is negative. 

The problem now is how to define and estimate the number N s of single NN 
collisions. Ideally, we would like them to be those independent collisions where produced 
particles do not reinteract at all with the surrounding environment. In perfectly central 
collisions, they supposedly are single NN collisions where both nucleons undergo exactly 
one collision occurring at the edge of the overlap region. On the other hand, in extreme 
peripheral nuclear collisions, they should reduce to one NN collision. In all other cases, 
they are tightly related to the NN collisions occurring at the edge of the overlap region 
where only one nucleon from either nucleus is involved, but their number cannot be 
defined in a clearcut way without a full dynamical model of the collision. However, we 
can resort to a definition interpolating the perfectly central and the extreme peripheral 
case and relying on the Glauber model. Such a definition might be: 

N s = mm[N 1(a) ,N m ] (3) 

where Nu a ) (Num) is the number of nucleons colliding once according to the Glauber 
model. We estimated the thus defined N s as a function of centrality by means of a 
Glauber Monte-Carlo calculation. The resulting N s /Ny/ ratio can be reasonably fitted 
for N w > 10 by (see fig. H]): 

N - 1 0.58 (4) 



N w AT™ 8 

Plugging (j3j) into ([2]) we have an expression of the normalized yield as a function 
of the number of wounded nucleons, that is centrality, depending on one unknown 
parameter A. This can be determined by matching the model to the measured value in 
the most central bin and then the centrality evolution is completely determined. The 
obtained curve is in impressive agreement with the data, as shown in fig. [5j the formula 
matches the experimental points to a high degree of accuracy. This is a clear evidence 
that the envisaged core-corona superposition is able to account for the strangeness 
undersaturation phenomenon. 

The same exercise can be repeated for open strange particles, the result being shown 
in fig. [5j It can be seen that the curves match the data in the most central bins, while 
they overestimate the measured points in most peripheral bins: this is likely due to the 
canonical suppression effect in the core which is not taken into account in the formula 
(J2J) but should indeed show up for peripheral collisions. Finally, we observe that other 
definitions of N s are possible (e.g. (iV 1(a ) + N 1 ( b - ) )/2) but they lead to similar results. 

5. Discussion and conclusions 



If our interpretation of production as a function of centrality is correct, several 
remarkable consequences are implied. First of all, the enhancement of relative strange 
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Figure 5. <\> (left panel) and hyperons (right panel) rapidity density per wounded 
nucleon as a function of participants normalized to pp collisions. Data points from 
STAR [T7|; solid lines are the predictions from core-corona superposition (see text). 



particle production going from peripheral to central collisions is mainly due to a 
geometrical effect of core-corona superposition. Canonical suppression plays a role only 
in the most peripheral collisions and it is possible that the SCV simply coincides with 
the core volume, what would be a nice fact itself. Secondly, the data supports evidence 
for a completely equilibrated hadron gas in the core throughout all centralities at RHIC, 
whose temperature is constant and equal to 165 MeV. 

The same conclusion is likely to apply to SPS too. The fact that there 75 ~ 0.85 
in central collisions [6], significantly lower than at RHIC, is related to the lower weight 
of the core compared to the corona. Indeed, as energy decreases, so does the freeze-out 
volume of the core and the multiplicity of particles stemming from it, while the the 
number of single NN collisions decreases only slightly, the NN cross section being slowly 
varying. This would nicely explain the mild increase of 75 as a function of centre-of- 
mass energy (see fig. E]), nevertheless a complete reanalysis of the data is compelling. As 
has been mentioned, early analyses of central collisions at top SPS energy based on this 
picture were fairly succesful [8], but peripheral collisions are indespensable to confirm 
this idea. In this respect, NA49 is going to update preliminary measurements [18] which 
were used to determine 75 in peripheral Pb-Pb collisions [19] . 

The constancy of T as a function of centrality (see fig. [3]) which was first observed 
by the STAR collaboration [20] is confirmed in our analysis of RHIC data [7] to a high 
degree of accuracy. This stunning independence of centrality is hard to reconcile with 
collisional thermalization, as pointed out in ref. [21] as it would require a dramatic 
dependence of hadronic reaction rates on temperature. Also, the Q yield is very difficult 
to reproduce in such an approach [22] unless invoking the existing of massive resonant 
degrees of freedom [23]. 
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